| Pembroke College 
Cambridge 


UNIVERSITY OF CAMBRIDGE 
PEMBROKE STREET, CAMBRIDGE 


ARTICLE 


Bubble Formation at an Orifice in 
a Viscous and in an Inviscid 


Liquid 
Author: Author: 
John Frank DAVIDSON B. O. G. SCHULER, 


1926 — 2019 


Transactions of the Institute of Chemical Engineers, 


Volume 38, pp. 144-154 & pp. 335-342 (1960) 


Contents 


Notation 


Bubble Formation at an Orifice in a Viscous Liquid 


. Introduction 
dd, Erteratyre survey: dues nee a a e a 
1:22, Present wörk ass 4- 2 rd rad ee bee e a al 


Experimental Methods 


Discussion and Results 


3.1. Bubble formation at constant flow rates... . 2.2222... 
3.1.1. SMALL FLOW RATES of ace deny do aes Be gael ae as 
3.1.2. LARGE FLOW RATES ero apar doe aeg 

3.2. Bubble formation with a constant pressure supply ......... 
3.2.1. SMALL FLOW RATES sara 2 a a eked 2. dl 
3.2.2. LARGE FLOW RATES. 2.0.4.6 4 2.446 45544 Rb 5 


Conclusions 


Acknowledgments 


Bubble Formation at an Orifice in an Inviscid Liquid 
Introduction 


Experimental Methods 


7.1. Apparatus for bubble formation ..................0. 


Constant Flow Regime 


8.1. Theoretical Considerations ................ 0008048 
8.2. Experimental Results . . 2.3222 a eee Ci eee hei 


Contents 


9. 


Constant Pressure Regime 

9.1. Theoretical Considerations a e A SS 
9.1.1. k SMALL; APPROACH TO CONSTANT FLOW ...... 
9.1.2. P = 0; THE CRITICAL FLOW ad ar daa 

9.2. Experimental Results 2 . . iaa & ra e a wee 


10. Conclusions 


11.Acknowledgment 


A. 


B 
C. 
D 


m 


Effect of the Perspex ring on bubble frequency 


. The effect of the gas momentum on bubble size 


Pressure necessary to generate a bubble in a viscous liquid 


. Pressure due to kinetic energy imparted to the liquid 


. Note on the formation of bubbles in an inviscid liquid 


John Frank Davidson. Obituary Notes 


List of Figures 


2 
2.2; 


3.1. 
3.2. 
3.3. 
3.4. 


3.0. 


8.1. 


8.2. 


9.1. 


9:2: 


9.3. 
9.4. 


9.5. 
9.6. 
9.7. 
9.8. 


9.9. 


A. 
B.1. 


Schematic diagram of apparatus used for formation of bubbles . . . 16 
E A be BT ta A Es Se Ee Boe a ee Dean 17 
Idealised sequence of events during bubble formation in a liquid .. 20 
Ciné photographs of bubble formation in a viscous liquid ..... . 21 
Comparison of theoretical and experimental bubble volumes for 
small constant gas flow Tales: Ltda e A A 23 
Formation of a bubble with coalescence of a small secondary bubble 
immediately after detachment . ... 2:2. 22. 2m rennen 24 
Assumed position of the bubble at detachment, using large orifice 
diametets: O A ze ne Pa A aah 25 
Experimental and theoretical bubble volumes for small constant flow 


TALES: er a re ie Meng Be Sa 45 
Experimental and theoretical bubble volumes for large constant flow 
PALES ge ee ee aa ee ner 4 47 
Bubble formation in groups of two with a constant flow rate of 
Tomales A A era AOE OTe Ba ed ie 49 
The theoretical relations between dimensionless bubble volume and 
PICSSUPCs << a o a ates Dae eS Ae u Wan 50 


The theoretical relations between dimensionless flow rate and pressure 51 


The theoretical relations between dimensionless bubble volume and 
HOw ahs Dal a high ag tec ak Cee ns Ae en ba ee cae eB a Sa 52 
Bubble volume as a function of pressure with air and water... . . 54 
Mean flow rate as a function of pressure, with air and water .... 55 
Bubble volume as a function of mean flow for constant pressure with 
is c.f erie AA gee is oS adios ee Peter le Bus dled 56 
Bubble formation at constant pressure with a mean flow rate below 
the critical value Re Sten ee A Sg wae Sark BR! A gs Re 57 
Bubble formation at constant pressure with a mean flow rate above 
thëcriticalk value Arte ad oles el Rew ein ge Ee ke Mes Be le 58 
Bubble frequency as a function of orifice ring diameter ....... 65 
Bubble sizes compared for carbon dioxide and air .......... 67 


List of Figures 


F.1. JOHN FRANK DAVIDSON, bromide print, 1993 


List of Tables 


3.1. 


3.2. 


3.3. 


3.4. 


9.1. 


Formation of Air Bubbles in Aqueous Glycerol at Large Constant 
Flow: Rates: ze oh nae Seles aoe ad Puls Gh ee tae Se a a vase 
Formation of Air Bubbles in Aqueous Glycerol at Small Flow Rates 
and Constant Pressure . . dl bad 2 La se ee oS 
Formation of Air Bubbles in Aqueous Glycerol at Large Flow Rates 
and Constant Pressure a 4035 a Fra ee 
Theoretical Calculation of Bubble Volume with constant Gas Pres- 
sure and Viscous Liquids. Results from ”Edsac 2”... ........ 


Formation of Air Bubbles in Water with Constant Gas Pressure . . 


Notation 


yr Pueeo ree 


2 . rv(3/ An) 3 

effective orifice area 
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distance of centre of sphere from wall 
3A/G 
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(constant) gas flow rate 


mean flow rate in "constant pressure” experiment 
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acceleration due to gravity 

depth of liquid seal 


orifice constant = flow rate 
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(pressure difference) 2 


parameter in equation (x) 
mass of liquid displaced by sphere 
mass of sphere 


Pı — pgh 
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-I Dt 
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mean pressure due to kinetic energy 


pressure in drum or pressure in the vessel below 
the orifice respectively 


pressure at radius R, 

Reynolds number 

radius of shell enveloping source of gas supply 
radius of bubble at any instant 

orifice radius 


distance of centre of bubble from point of gas 
supply, vertical distance moved by the bubble 


time, time bubble has been growing 
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time at detachment 
radial velocity of liquid surrounding bubble 
bubble volume at any instant /time t 
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Part I. 


Bubble Formation at an Orifice in a 
Viscous Liquid 


Abstract 


A theory of bubble formation based on the motion of a bubble in a viscous liquid 
has been developed. The theory gives the volume of gas bubbles formed at an 
orifice in a viscous liquid for both constant gas flow and constant gas pressure. 
Experiments were carried out with liquids of high viscosity (500 cP to 1040cP)?. 
Good agreement with theory was obtained over a large range of gas flow rates 
(OmLs~! to 50mLs”!). 


‘Equals in SI units 500 mPas to 1040 mPas, which equals approximately the viscosity of must 
and fuel oil respectively 


1. Introduction 


1.1. Literature survey 


When air is blown steadily through an orifice submerged in a liquid it emerges as 
bubbles which are formed periodically over a wide range of flow rates. This paper 
describes measurements of the frequency of the bubbles and also the calculation of 
the frequency from first principles. 

A comprehensive literature survey on the subject of the formation of bubbles at 
submerged orifices was described by HUGHES, HANDLOS, EVANS, and MAYCOCK, 
[15], in 1953. Several publications dealing with various aspects of the problem have 
subsequently appeared. DAVIDSON and AMICK, [7], studied the effect of the volume 
of the orifice-chamber on the size of the bubbles, whilst HELSBY and Tuson, [14], 
photographed different modes of bubble formation. BENZING and MYERS, [1] 
formed bubbles of air in water at very small gas flow rates, and found that the 
diameter of the bubble could conveniently be calculated by means of a correlation 
which did not include a viscosity term. QUIGLEY, JOHNSON, and HARRIS, [23] 
varied the viscosity of the liquid between 1 and 400 cP and the density between 1.00 
and 1.57 g mL”*, and concluded that the effect of both these variables on bubble-size 
was small. Similarly the effect of surface tension was negligible at higher flow rates. 
The controlling factors appeared to be orifice diameter and gas flow rate. LEIBSON, 
HOLCOMB, CACOSO, and JACMIC, [16,17], used larger flow rates and found that 
at high Reynolds numbers the formation was not periodic: the size distribution 
fitted a logarithmic normal probability distribution. CALDERBANK, [3,4], used 
orifices and slots ranging in diameter from = in (0.16cm to 0.64cm), and has 
suggested that for air bubbling into water at flow rates between 20 and 250 mL st, 
the frequency is nearly constant within the range 15-20 bubbles s”!. 

Further relevant data have been published in a series of papers by SIEMES and his co- 
workers, [24-29]. Of particular interest is the work of SIEMES and KAUFMANN, [28], 
who studied the periodic formation of bubbles at single submerged orifices under 
various physical conditions. They arrived at important conclusions: 


(i) With small gas flow rates and with liquids of low viscosity the difference 
in volume between the bubble formed at any particular flow rate and that 
formed at infinitely small flow rates using the same orifice is independent 
of the viscosity, surface tension, and density of the liquid. The controlling 
factors are gas flow rate and orifice diameter. 
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(11) For liquids of high viscosity the effect of surface tension and density is still 
negligible, but the viscosity becomes an additional factor controlling bubble 
size. It is at this point worth noting that for viscous liquids the volume of the 
bubble formed at infinitely small gas flow rates is negligibly small compared 
with the total volume formed at higher flow rates. 


SIEMES and KAUFMANN suggested the following mechanism of bubble formation. 
They postulated that a bubble forms in two distinct stages. For inviscid liquids the 
initial stage is presumed to be exactly analogous to the formation of a bubble when 
the gas flow rate approaches zero. The second stage is presumed to begin when the 
buoyancy forces exactly balance the surface tension forces and the bubble proceeds 
to detach itself. The detachment time is presumed to be constant, irrespective 
of the gas flow rate. Gas flows into the bubble during this time and hence the 
influence of gas flow rate on bubble volume can be derived. For viscous liquids a 
similar mechanism was suggested. Initially the bubble is supposed to grow to a 
definite volume, which in this case is to some extent dependent on the viscosity of 
the liquid but, once again, not on the gas flow rate. As before, the bubble then 
requires a constant period to detach itself, so that the flow rate affects the bubble 
volume. The actual time for the process of detachment to occur was calculated to 
be several hundredths of a second. 

The experimental conditions described by SIEMES and KAUFMANN suggest that 
their apparatus delivered a constant flow of gas, not as in the case where gas flows 
from a chamber at constant pressure through an orifice into the forming bubble. 
The difference between the two cases will be considered in greater detail below. 
Finally it should be noted that the mechanism put forward by SIEMES and KAUF- 
MANN is based entirely upon a description of a suggested physical process. The 
explanation offers no way of predicting the time of detachment, and for this reason 
the theory cannot be regarded as being derived from first principles. 


1.2. Present work 


The present paper deals with the formation of bubbles in liquids of high viscosity 
and with gas flow rates between 0OmLs~! to 50mLs~!. Two distinct situations are 
considered, namely:— 


1. In the simplest case there is a constant flow of gas throughout bubble formation 
as when the gas is being fed through a long thin capillary. The surface of 
the bubble, as it grows, is assumed to be spherical, so that the radius of the 
bubble at any instant of time, t, is known. As the bubble grows it tends to 
move upwards, the force being due to buoyancy. This force has to overcome 
two resistances, namely, viscous drag of the STOKES kind!, and liquid inertia 


ISTOKES's law is, although it can be derived theoretically, an empirical law for the frictional 
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due to the fact that as the bubble accelerates upward some of the liquid also 
accelerates. The balance between the upward force and the two resistances 
leads to a differential equation governing the motion, and the solution gives 
the distance of the centre of the bubble above the orifice at any time t. 
Initially, when the bubble is small, its centre is at the orifice. The formation 
is complete when the bottom of the bubble reaches the orifice, and at this 
instant the bubble detaches and the process then repeats itself. 


2. In the second case the bubbles are formed above a hole in a horizontal plate, 
with liquid above to a finite depth and gas below at a constant pressure. In 
this case, which is of greater practical importance, the gas flow through the 
orifice varies during bubble formation. This has the effect of making more 
complex the equations governing the motion of the bubble. The basis on 
which the equations are formulated is, however, similar to that described in 
section 1.1 above. 


The theory bears out several of the findings of SIEMES and KAUFMANN. For a 
constant flow system the bubble size depends upon flow rate and viscosity, whilst 
surface tension has no effect. However, with regard to the effects of liquid density 
and orifice diameter, the theoretical results are in conflict with the conclusions 
Of SIEMES and KAUFMANN. The liquid density is shown to affect the size of 
the bubble whilst the influence of the orifice dimensions is negligible. The orifice 
diameter becomes important only when the effective area is small and the gas 
flow rate relatively large. Under these conditions the gas receives a considerable 
amount of momentum which affects the motion of the bubble and hence its final 
volume. For a constant pressure system, both the surface tension of the liquid and 
the orifice discharge coefficient have to be included in the equations governing the 
motion of the bubble. The predicted values of bubble volume for both constant flow 
and constant pressure systems were found to be in good agreement with results 
obtained experimentally. 


force — also called drag force — exerted on spherical objects with very small Reynolds numbers 
in a viscous fluid. It was derived by GEORGE GABRIEL STOKES in 1851 by solving the Stokes 
flow limit for small Reynolds numbers of the NAVIER-STOKES equations. 

The force of viscosity on a small sphere moving through a viscous fluid is given by: 


Fa = 6rpav 


where: 

Fa is the frictional force — known as STOKES’ drag — acting on the interface between the fluid 
and the particle, u is the dynamic viscosity of the fluid, a is the radius of the spherical object, 
and v is the flow velocity relative to the object. 
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2. Experimental Methods 


A schematic diagram of the apparatus is shown in Figure 2.1. Gas from either a 
compressor or a cylinder was passed into a 25 L vessel, A, which acted as a buffer. 
The pressure in A was indicated by the Manometer C and could be maintained at 
any desired value by manipulating valve B. In most cases the pressure was kept 
at about 1.5 atmospheres (1.52 bar) absolute. Next the gas passed through an 
EDWARDS throttle valve, D, by means of which the flow rate could be controlled. 
The threeway tap, E, enabled the gas to pass either directly to the drying tube, 
H, and the Rotameter, I, or first through a soap bubble meter, F. J was a 45 L 
drum from which the gas passed through the orifice, K, into the liquid contained 
in a Perspex cylinder of 14.7cm diameter. The latter was attached directly to the 
drum and was surrounded by a thermostat. L and M were helical coils immersed in 
the same liquid as was circulated through the thermostat surrounding the Perspex 
cylinder. Gas temperatures were measured at N and O. These were always within 
0.2°C of the temperature of the liquid in which the bubbles were forming. The 
temperature of the liquid itself did not vary by more than 0.1°C during a run. The 
frequency of bubble formation was determined by stroboscopic illumination and the 
pressure in J was given approximately by the water Manometer, G. More accurate 
pressure determinations were made by means of a tilting Manometer connected 
to the drum J. A FERRANTI Viscometer was used to determine the viscosity of 
the liquid, and the surface tension was measured by the drop weight method. The 
liquids used in most of the experiments were aqueous solutions of glycerol for which 
density values were obtained from the International Critical Tables. In the few 
cases where viscous hydrocarbon oils were used, the densities were determined with 
the aid of a Pyknometer. Provided the orifice was immersed in liquid to a depth of 
more than a few bubble diameters, this depth had no effect upon bubble frequency 
but merely altered the pressure necessary to generate the bubbles. The shapes of 
the orifices used in the experiments are shown in Figure 2.2. 

Figure 2.2(a) represents the type of orifice used in the constant pressure systems. 
The nozzle was made of brass and had a Perspex ring fitted as shown. The purpose 
of the ring was to eliminate as far as possible the effects of the circulation of the 
liquid (see Appendix A). In these experiments the pressure drop across the orifice 
was of the order of a few centimetres of water. 

For the experiments in which a constant flow throughout bubble formation is 
required, it is necessary to have a large pressure drop across the orifice, which can 
be achieved simply by making the diameter of the orifice very small. This, however, 
has the disadvantage that the neck of the bubble is short at detachment and the 
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Figure 2.1.: Schematic diagram of the apparatus used for the formation of bubbles 


bubble subsequently is in such a position that for certain critical flow rates the next 
bubble, which rapidly expands in the initial stages, coalesces with it. The larger 
the orifice diameter, the greater the gas flow rate at which coalescence first occurs. 
For large flow rates, a small orifice has the additional disadvantage of imparting a 
large momentum to the gas, which, although it can be approximately accounted 
for, complicates the theoretical calculations. On the other hand the theoretical 
boundary conditions are more closely obeyed if the neck of the bubble is short. 
Bearing these considerations in mind three types of orifice were eventually used: 


(i) 


(i) 


(iii) 


For very small flow rates (0mLs~! to 2mLs”!) the orifices were essentially 
the same as those used in the constant pressure experiments, except that the 
orifice diameters were small and hence the pressure drops large. Pressure 
fluctuations within the bubbles due to changes in radius of curvature had a 
negligible effect on the flow, especially as the latter depends on the square 
root of the pressure difference. 


Experiments were also performed in which the drum J was by-passed and 
the gas fed directly into the liquid through a length of capillary. A ring, as 
shown in Figure 2.2(c), was attached to the end of the capillary and, provided 
the diameter of the capillary was small, the results were the same as those 
obtained when the small diameter orifices, described above, were used. 


For larger constant flow rates the orifices used were of the type shown in 
Figure 2.2(b). A piece of sintered brass was welded to the lower side of the 
orifice which in other respects was similar to that shown in Figure 2.2(a). 
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Perspex ring 
Brass nozzle 


Capillary 
(a) Rubber washer 
Sintered brass 
plate 
u Plate attached 
to drum «) 


JAV. 


(b) 


Figure 2.2.: Orifices used 


(a) equals constant pressure systems 
(b) equals constant flow systems (large flow rates) 
(c) equals constant flow systems (small flow rates) 


The diameters of the orifices were of the same order as those used in the 
constant pressure experiments. A high pressure drop could be maintained 
across the sintered brass, resulting in a constant flow of gas into the bubble. 
The volume of the cylindrical cavity between the sintered brass and the 
upper plane of the nozzle was small compared with the volume of the final 
bubble and its effect was neglected. The purpose of the relatively large orifice 
diameter was two-fold. Firstly the bubble during detachment had a slightly 
longer neck. This meant that when the bubble detached, it was far enough 
away from the orifice to avoid contact with the next bubble. In this way 
the condition assumed in the theory, namely the periodic formation of single 
bubbles, was satisfied. Secondly the large cross-sectional area of the orifice 
kept the momentum of the gas low. 


To help elucidate the mechanism of the formation of a bubble it was decided to 
take ciné photographs showing the entire formation. The design of the vessel in 
which the bubbles were being formed was slightly modified in order to eliminate 
distortion of the picture due to the curvature of the Perspex cylinder, and the 
stroboscope frequency adjusted until one bubble appeared to be forming every 15s 
to 20s. A photograph of the bubble was then taken every 0.75s using an externally 
triggered 16mm ” Pathé” camera. ” KODAK Plus X” proved to be a convenient film 
for an exposure of 0.2s. The frequency of the stroboscopic illumination was usually 
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between 9 and 21 flashes st, so that during each exposure between 1 and 5 pictures 
were superimposed upon each other. The light output of the DAWE Strobotorch 
used was approximately constant, irrespective of the frequency, so that no change 
of film was necessary when changing from one frequency to another. 
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3. Discussion and Results 


The following theory attempts to describe the way in which gas bubbles are formed 
at an orifice in a liquid. It is assumed that the gas is supplied at a point source 
within the liquid and that the bubbles, as they form, are spherical. Each bubble 
is assumed to begin with its centre at the point source and the upward motion 
is determined by a balance between the upward force due to buoyancy and the 
drag forces due to viscosity and inertia. Now initially the upward velocity of the 
centre is small since it starts from rest, and because the bubble is expanding the 
underside has a downward velocity and the bubble therefore continues to envelop 
the source of gas. Subsequently, when the bubble is larger, its base is brought to 
rest when the outward velocity of the bubble surface relative to its centre is equal 
to the upward velocity of the centre due to buoyancy. Thereafter the base has a 
net upward velocity and detachment takes place when it reaches the source of the 
gas. The idealised sequence of events is indicated in Figure 3.1, and this picture of 
formation and detachment was the basis of all the theoretical calculations, although 
in some cases it was assumed that detachment was delayed until the base of the 
bubble was somewhat above the source of gas. 

The theory given below is in two sections. The first section deals with the 
formation of bubbles at a constant gas flow rate and the second with the case 
where gas flows through an orifice from a vessel at constant pressure. Each section 
is subdivided to distinguish between small and large flow rates. The assumptions 
made and the relevant results obtained both by calculation and experiment are 
given in each case. Obviously the exact conditions assumed by the theory cannot be 
realised experimentally, especially if the bubbles are formed just above horizontal 
plates. Photographs (Figure 3.2) show that the bubbles are pear shaped initially, 
this being due to the fact that the base of the bubble cannot move downwards. 
Also, some of the boundary conditions, to be discussed in detail below, are not 
exactly obeyed. Nevertheless, the predicted values are in good agreement with the 
experimental results. 
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(a) (b) (e) 


Figure 3.1.: Idealised sequence of events during bubble formation in a liquid. Gas 
supplied from a point source P 


(a) = underside of bubble moving down 
(b) = bottom point of bubble at rest in lowest position 
(c) = bottom point of bubble reaches gas source 


3.1. Bubble formation at constant flow rates 
3.1.1. SMALL FLOW RATES 
The theory is based on the following assumptions: 
1. The bubble is spherical throughout formation. 
2. Circulation of the liquid is negligible, so that the liquid surrounding the 
orifice is at rest or stagnant respectively when the bubble starts to form (see 


Appendix A). 


3. The motion of the bubble is not affected by the presence of another bubble 
immediately above it. 


4. The momentum of the gas is negligible (see Appendix B). 


5. The bubble is at all instants moving at the STOKES velocity appropriate to 
its size. 


Consider the motion of a bubble forming at a point source in an infinite liquid 
under the above conditions. In view of assumption 5., the velocity of the centre of 
the bubble, v, at time t after the start will be given by the equation: 


(3.1) 
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Figure 3.2.: Ciné photographs of bubble formation in a viscous liquid 


Gas flow rate = 24mLs”! 
Viscosity = 711cP (711 mPas) 
Density = 1.25gmL7! 

Radius of orifice = 0.096 cm 
Bubble volume = 2.5 mL 


Also, if G is the gas flow rate, then, V, the bubble volume is given by 


4 3 
G=V= > (3.2) 
where r is the bubble radius. Hence 
2 
2g (3Gt\? ds 
en, Bat 3.3 
“=p ( An ) dt oy 


where s is the distance between the centre of the bubble and the point of gas supply. 
Integrating Equation 3.3 we have 


_ 2g 3G os 
“A Xd 


wI 


(3.4) 
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The bubble will detach when s = r, and hence at detachment, eliminating t from 
Equation 3.2 and Equation 3.4, 


V= a ES . (3.5) 


Equation 3.5 gives the volume of the bubble in terms of the properties of the liquid 
and the gas flow rate. A comparison of theoretical and experimental results for 
aqueous solutions of glycerol is presented in Figure 3.3. Runs were also done for 
similar flow rates using viscous hydrocarbon oils of varying surface tension. In each 
case Equation 3.5 was closely obeyed. A change in surface tension alone thus has 
no effect on the bubble size. 

At a gas flow rate of about 1.5mLs”!, the apparent break in the sequence of 
points corresponding to a viscosity of 515 cP (515mPas) is due to the coalescence 
of a small secondary bubble with the main one as the latter leaves the orifice 
(Figure 3.4). The discontinuities in the curves of SIEMES and KAUFMANN are 
probably due to this phenomenon. 

The assumption that the bubble is at all times moving at its STOKES velocity is 
at best only an approximation. GARNER and HAMMERTON, [10], measured the 
ultimate velocities of small bubbles and came to the conclusion that the drag on 
the bubble depends on the Reynolds number, defined as Re = ain where dp is the 
diameter of the equivalent sphere. For glycerol, STOKES law was found to apply at 
Re < 0.002. At higher values (0.023 < Re < 0.2) there was agreement with the 
HADAMARD correction, while at still higher values a transition took place in the 
opposite direction, the drag becoming considerably larger than the STOKES drag 
at Re = 10, [20]. At these higher Reynolds numbers, the shape of the bubble tends 
to deviate from the spherical form. 

If a bubble is formed as described above, the Reynolds number will change with 
time, the range for these small flow rates being from zero to about 1.5, and the 
STOKES drag was in the first instance chosen because it represents an approximate 
average. 

It should also be noted that the experimental results deviate to some extent from 
the theoretical predictions as the flow rates are increased. This is probably due to 
the fact that the inertia of the surrounding liquid was ignored. 


3.1.2. LARGE FLOW RATES 


As mentioned above, for large flow rates it cannot be assumed that the inertia 
terms are negligible, and in the theory which follows, an attempt has been made 
to allow for the effect of inertia. 

If a sphere moves in a direction perpendicular to a wall the virtual mass, [19], of 


the sphere is given by M, + M A m 
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BUBBLE VOLUME (mi) 


o os ro es 29 25 


GAS FLOW RATE (mi/») 
Orifice radius — 0-0334 cm + = 764 cp 
Ring radius = 1:10 cm ) = 695 cp 
Liquid seal = 5-0 cm A = 626 cp 
O = 1040 cp x = 515 cp 
V = 915 cp 


Figure 3.3.: Comparison of theoretical and experimental bubble volumes for small 
constant gas flow rates. The lines give the theoretical values corresponding to the 
nearest set of experimental points 
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Figure 3.4.: Formation of a bubble with coalescence of a small secondary bubble 
immediately after detachment 


If the sphere is a bubble, M, is negligible, and the virtual mass is M - (3 + ia). 
When a = c, the virtual mass is uM and this was taken as an average value of the 
virtual mass for a bubble forming above a horizontal plate. 
It should here be mentioned that the above considerations regarding the virtual 
mass apply only if the liquid is completely inviscid. This condition is not satisfied 
in the present work and the calculations below can only be regarded as giving 
the order of magnitude of the effect due to inertia. However, the viscous forces 
outweigh the inertia forces. Assumptions 1. to 4. listed in subsection 3.1.1 remain 
unchanged and the equation governing the motion of the bubble becomes 

11 d(Mv) 


Oa a + 6rvpro, (3.6) 


where p is the density of the liquid. Substitution in Equation 3.6 of pGt for M and 
1 

3Gt 3 : 

Er for r gives 

3Gt5 


11 
+ 6 + 6av r -v. (3.7) 


dv 


dt 


11 
ea 
g 16 


Equation 3.7 applies if the bubble is assumed to detach when s = r, as when the 
bubble is forming at a point source. If, however, the radius of the orifice is large, 
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Figure 3.5.: Assumed position of the bubble at detachment, using large orifice 


diameters 


then a better value of s at detachment is s = r + ro, where ro is the radius of 
the orifice. At detachment the bubble would then be in the position indicated in 
Figure 3.5, and Equation 3.7 has to be modified to allow for the initial volume of 


the bubble, Vo = sty at t = 0. The modified equation is 


11 dv 11 33 1 
g- (GE+ Vo) = 15: (GEH Vo) + eG + bry > (Gt+ Vo)? 


dt An 


Rearranging and substituting 


and V for Gt + Vo, Equation 3.8 becomes 


Qe ome q 
dt Vo yl T 


(3.8) 


(3.9) 
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which, on integrating from V = Vo to V and assuming that at V = Vo, v = 0, gives 


ds 489 
v= a A, 
dV 11G 
vs 5V3 20 BV 120V-3  120V” L (3.10) 
D- D DP. De p Dé V.epv w 
where 
_ 489 pył 
11G 
5 4 1 
Vè 5V _ 20V 60V3 120V 120 
D D2" D Dt ` Dd D6 
and A 
D=— 
G 
A further integration from V = Vo, where s = 0, to V gives 
_ 48g 
Te 
5 V 
3V3 15V5 20V 90V: 360V3  120In(V) 
5D 4D? D Dt ` D’ Dp | 
0 
n= V 
: sen 
A A 3.11 
= |m) +3 2, a (3.11) 
= E 


Equation 3.11 enables a graph of s versus t to be constructed and the time for the 
bubble to reach a volume such that s = r + ro can thus be determined. This time 
multiplied by the gas flow rate then gives the volume of the bubble that detaches 
and moves away. Theoretical and experimental values of the volumes of bubbles 
are compared in Table 3.1. 
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Flow rate Kinematic Bubble Bubble Theory 
(mL s~+) viscosity volume volume Expt. 
(Stokes) (experiment) (theory) 
(mL) (mL) 
18 7.9 2.1 1.9 0.90 
24 7.9 2.6 2.5 0.96 
33 7.9 3.2 3.3 1.03 
15 6.5 1.6 1.6 1.00 
26 6.5 2.4 2.4 1.00 
35 6.5 3.0 3.1 1.03 
14 5.4 1.3 1.3 1.00 
20 5.4 1.8 1.7 0.95 
26 5.4 2.2 2.2 1.00 


Table 3.1.: Formation of Air Bubbles in Aqueous Glycerol at Large Constant 
Flow Rates, rọ = 0.096 cm 


3.2. Bubble formation with a constant pressure 
supply 


3.2.1. SMALL FLOW RATES 


For each orifice considered, the flow of gas is proportional to the square root of the 
pressure difference across the orifice. 

Suppose now that an orifice of this kind is submerged in a liquid and that a 
bubble is forming on the upper surface. Let the gauge pressure on the lower side 
of the orifice be P,. The gauge pressure on the upper side of the orifice will be 
pgh + = — pgs, where h is the head of the liquid and o the surface tension. 

The rate at which gas will flow into the forming bubble will thus at any instant be 
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given by: 
20 dV 
G=k-4/ P, — pgh — ^ => 3.12 
yP. pgh — — + pgs = 7 > (3.12) 
where k is an experimental constant depending on the orifice and the gas flowing 
through it. 


The effect of the kinetic energy of the liquid on the pressure within the bubble has 
been neglected; an expression derived in Appendix B shows that this effect is small. 
It is also shown in Appendix C that the viscosity of the liquid has only a small 
effect upon the pressure within the bubble. 

For small gas flow rates, r and s will be small and as a first approximation the 
term pgs can be neglected and Equation 3.12 becomes 


2 dV 
=k. —poh— 2 = E. (3.13) 
Substituting r = (3V/Ar)3 and writing P for Pı — pgh, Equation 3.13 becomes 


dv dn \ 3 
ei Pr A 
dt e j (= 


which can be integrated by substituting 


Ar 3 20 
Del — 
ASE 3V P 


3 
Arro 
3 


and using the boundary condition V = 


Bet (2) 


at t = 0 to give 


T 
3 P 


5 1 
. fiano . (soc 0 + E sec? 0 + a - Sec 0) (3.14) 


+ -In tan (3 + 3) \ (3.15) 


the limits of integration being 


0 = l - 
fa 
and 
1 
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V and r can thus be calculated for any value of t. 

The assumptions listed in subsection 3.1.1 are applicable, and Equation 3.1 and 
Equation 3.14 enable a plot of v versus t to be constructed, from which by means 
of graphical integration the distance s can be computed for any time, t. A plot of 
s versus r+ ro then gives the volume of the bubble when s = r + ro. A comparison 
of experimental and theoretical results for air and aqueous solutions of glycerol is 
given in Table 3.2. 

The orifice coefficient k was determined with air on both sides of the orifice plate; 
k was assumed to be the same with air bubbling into liquid. 

Table 3.2 shows that the theoretical and experimental values of bubble volumes are 
in reasonable agreement in view of the assumptions made in deriving the theory. In 
these constant pressure experiments, the flow rate is a dependent variable, and it 
is therefore appropriate to compare theoretical and experimental values of the flow 
rate. This comparison is given in Table 3.2, and the last column of figures shows 
good agreement between theory and experiment, except for the first figure, 1.66. 
In this case the pressure, 4860 g cm”! s? (4.86 x 1073 mbar), is almost exactly equal 
to a and it is thought that in the experiments there was an appreciable delay 
period between the completion of the bubble and the beginning of the expansion 
of the next bubble. In the theory it was assumed that there was no delay period 
between bubbles, so the theoretical flow should be higher, as shown in the first 
line of Table 3.2. This leads to the concept of a critical flow for continuous bubble 
formation which will be discussed in a subsequent paper. 

Experiments were also carried out using carbon dioxide instead of air. The orifice 
coefficient was redetermined, and as in the case of air, the size of the bubbles 
formed was in good agreement with theory. 
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Bubble volume Flow rate 
— A Á A ==, ee A A, 
Orifice Orifice Kinematic Surface Pressure Expt. Theory Theory Expt. Theory Theory 

constant k radius viscosity Tension  (gcm”?!s?) (mL) (mL) Expt. (mL s~!) (mL s~') Expt. 
(mL em2g 2) cm (Stokes) (dyn cm7!) 

0.0414 0.0260 7.8 63.9 4860 0.46 0.39 0.85 1.7 2.82 1.66 

0.0414 0.0260 7.8 63.9 7080 0.51 0.49 0.96 3.1 3.86 1.24 

0.0414 0.0260 7.8 63.9 8930 0.55 0.53 0.96 3.6 4.24 1.18 

0.0414 0.0260 3.6 64.7 6110 0.26 0.24 0.92 2T 3.12 1.15 

0.0414 0.0260 3.6 64.7 8290 0.30 0.28 0.93 3.4 3.78 1.11 

0.0901 0.0323 4.6 64.5 4530 0.56 0.47 0.84 5.1 5.88 1.15 

0.0901 0.0323 4.6 64.5 6920 0.67 0.57 0.85 Tel 7.60 1:07 


Table 3.2.: Formation of Air Bubbles in Aqueous Glycerol at Small Flow Rates and Constant Pressure 
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Bubble volume Flow rate 
—  áK — er “OO 
Kinematic Surface Pressure Expt. Theory Theory Expt. Theory Theory 

viscosity Tension (gem”?s?) (mL) (mL) Expt. (mL s”!) (mL s~) Expt. 
(Stokes) (dyn cm7!) 

7.4 63.9 2250 4.5 4.1 0.91 32 42.7 1.33 

7.4 63.9 3120 5.0 4.7 0.94 42 50.0 1.19 

7.4 63.9 5360 5.8 5.9 1.02 57 64.9 1.14 

6.2 64.0 1750 3.9 3.4 0.87 29 37.8 1.30 

6.2 64.0 2310 4.1 3.8 0.93 35 43.2 1.24 

6.2 64.0 3350 4.6 4.5 0.98 46 52.3 1.14 

5.4 64.1 2680 3.9 3.8 0.97 38 45.8 1.20 

5.4 64.1 3650 4.3 4.4 1.02 47 53.6 1.14 

5.4 64.1 4940 4.7 5.0 1.06 56 61.8 1.10 


Table 3.3.: Formation of Air Bubbles in Aqueous Glycerol at Large Flow Rates and Constant Pressure 


Orifice constant k = 0.856 mL cm? g73 
Orifice radius rọ = 0.0975 cm 
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3.2.2. LARGE FLOW RATES 


For large flow rates the assumption that the change in hydrostatic head does not 
affect the flow of gas into the bubble is no longer valid. The bubbles formed are 
large and the effect of the term pgs on the pressure difference becomes appreciable. 
Also, as in subsection 3.1.2, the inertia terms become important, and must be 
allowed for in the equation of motion of the bubble. 

The two simultaneous equations to be solved are Equation 3.12 and 


SA de. Tie ds. ds dv 
= rv - -— 4 | 3 
V, = 61 v ES dE (3.16) 


16 d ` 16 dt dt 
Equation 3.12 and Equation 3.16 cannot be solved analytically and values of s, V, 
and as were calculated on the electronic computer ” Edsac 2” for regularly increasing 
values of t, and Table 3.4 gives a typical set of results in detail, showing how s and 
V vary with time. The initial values of s, V, and ds for t = 0 were taken to be 0, 
Anré 


3 
was taken to be the total volume at s = r+ rg minus the initial volume ate = Vi. 


The volume Vp remained to form the nucleus of the next bubble. Table 3.3 gives a 
comparison of theoretical and experimental values of bubble volume, and of flow 
rate. The agreement between theory and experiment is satisfactory, except at the 
lowest pressures. 


, and 0 respectively. The volume of the bubble moving away from the orifice 
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t S V @ 

(s) (cm) (em) 2 

0.000 0.0000 0.0039 >1.000 
0.005 0.0051 0.1812 >1.000 
0.010 0.0181 0.3697 0.394 
0.015 0.0380 0.5617 0.240 
0.020 0.0642 0.7564 0.162 


0.025 0.0964 0.9537 0.126 
0.030 0.1341 1.1537 0.103 
0.035 0.1773 1.3565 0.086 
0.040 0.2257 1.5622 0.074 
0.045 0.2792 1.7711 0.064 


0.050 0.3376 1.9832 0.055 
0.055 0.4008 2.1988 0.049 
0.060 0.4688 2.4181 0.044 
0.065 0.5414 2.6411 0.038 
0.070 0.6187 2.8680 0.035 


0.075 0.7004 3.0990 0.033 
0.080 0.7866 3.3342 0.029 
0.085 0.8772 3.5737 0.027 
0.090 0.9722 3.8176 0.026 
0.095 1.0716 4.0662 0.023 


0.100 1.1752 4.3193 0.021 


Table 3.4.: Theoretical Calculation of Bubble Volume with constant Gas Pressure 
and Viscous Liquids. Results from ” Edsac 2” 


v = 7.4 St (u = 929 mPas or 929 cP) 

P = 2250 gcm! s7? (2.25 x 107° mbar) 
o = 63.9 dyn cm? (0.639 mN cm~?) 

k = 0.856 mL cm? g2 

ro = 0.0975 cm 

p = 1.256 g mL“! 
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4. Conclusions 


The volume of a bubble formed at an orifice submerged in a viscous liquid can be 
calculated by assuming that the orifice acts as a point source of gas supply and 
that the motion and position of the bubble is governed by viscous, inertia and 
buoyancy forces. 

Two cases can be distinguished, namely constant gas flow rate and constant supply 
pressure of the gas. In both cases the forces due to inertia are negligible if the 
overall gas flow rate is small. 

If the gas flow rate is constant, then the volume of the bubble that detaches is a 
function of the gas flow rate and kinematic viscosity of the liquid. The dimensions 
of the orifice are of minor importance. 

For constant gas pressure, however, the orifice dimensions become important and 
the volume of the bubble is a function of a constant depending on the size and 
shape of the orifice, the density of the gas, the kinematic viscosity and surface 
tension of the liquid and the pressure of the gas supply. In this case of constant 
pressure supply, the average gas flow rate is a dependent variable which is affected 
by all the factors which influence the bubble volume. 

The effects of the physical variables are summarised below. 


Viscosity. 
According to STOKES’ law the velocity of a bubble rising in a viscous liquid is 
inversely proportional to the viscosity and since the volume of a bubble forming 
above an orifice depends on the time it is in contact with the orifice, the viscosity 
has a major effect on bubble size. 


Gas flow rate. 
At any particular viscosity, the frequency of bubble formation is nearly constant, 
and changes in gas flow rate merely change the volumes of the individual bubbles. 


Orifice dimensions. 
The orifice dimensions are not important for constant gas flow rates. With constant 
pressure gas supply, however, the flow through the orifice is proportional to its 
cross-sectional area, making the latter very important. 


Liquid density. 
With constant gas flow rates, an increase in the density of the liquid has the effect 
of increasing the velocity with which the bubble rises and a smaller bubble results. 


4. Conclusions 


With constant pressure the density of the liquid also has an effect on the gas flow 
rate into the bubble. As the bubble rises, the hydrostatic head, which depends on 
the density of the liquid, decreases. The pressure in the bubble becomes less, and 
the gas flow rate into it increases. 


Surface tension. 
If the gas flow rate is constant, the surface tension has no effect other than 
that due to the small forces arising from contact round the edge of the orifice. 
With constant gas pressure, however, the surface tension has an appreciable effect 
on the pressure in the bubble and so to some extent governs the flow into the bubble. 


Gas density. 
Gas density is important only with constant pressure supplies, where it affects the 
discharge of gas through the orifice. The effect of the momentum of the emerging 
gas is small, over the range of variables considered in this paper. 
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Part Il. 


Bubble Formation at an Orifice in an 
Inviscid Liquid 


Abstract 


The periodic formation of bubbles due to the flow of gas into an inviscid liquid has 
been studied both experimentally and theoretically. Two systems were studied: 


(i) The simplest case is when gas flows through an orifice into the liquid at a 
constant volumetric rate, so that the movement and growth of the bubble do 
not influence the gas flow rate. In this case a simple theory, based on first 
principles, gives a rough estimate of the relation between bubble volume and 
flow rate. 


(ii) When gas passes into the forming bubble from an infinite vessel at constant 
pressure, the theory has to be modified to allow for the change in gas flow 
rate into the bubble as it forms. Nevertheless, the relation between gas flow 
rate, G, and bubble volume V is similar to the relation for case (i). 


However, when the pressure in the vessel is reduced, G and V are suddenly reduced 
to zero at a certain value of the vessel pressure, so that there are critical values of 
G and V. These critical values of pressure, flow, and bubble volume are in rough 
agreement with experimental values, and it is hoped that this theory will throw 
some light on the phenomenon of ”dumping” on sieve plates. 


6. Introduction 


This paper describes an experimental and theoretical investigation of the formation 
of gas bubbles at an orifice in an inviscid liquid. The main purpose of the paper is 
to describe a theory of bubble formation by means of which the size and frequency 
of bubbles can be calculated. This attempt has been only partly successful in 
that the size and frequency of the bubbles do not always agree with experiment, 
but the theory is nevertheless valuable in assessing the approximate effect of the 
various parameters. In particular, the theory is valuable in predicting that there is 
a critical gas flow rate for continuous bubble formation in the so-called ” constant 
pressure” experiments. In these experiments, air was passed through a single hole 
of radius rọ in a horizontal plate into a finite depth of liquid above the plate. It is 
clear that for bubbles to form, the pressure P, in the vessel below the plate must 
exceed a minimum value equal to the sum of the static head of liquid h above the 
plate plus a surface tension term 2, The theory shows — and experiments confirm 
— that if the vessel below the orifice plate is large, then as soon as P, is greater 
than pgh + a bubbles of a finite size will form continuously, so that the gas flow 
rate suddenly increases from zero to a finite amount. Here p is the liquid density 
and g the acceleration due to gravity. This finite gas flow is of the same order of 
magnitude as the critical gas flow per hole on a sieve plate, the so-called ” dumping 
rate”, [2,18]. 

The paper also deals with studies of bubble formation with a constant gas flow 
rate, the constant flow being obtained by feeding the air into the liquid through 
a restriction such as a capillary or sinter, so that there is a large pressure drop 
between the gas reservoir and the feed to the bubble. This arrangement is of less 
practical importance, but is simpler to treat theoretically and is therefore useful in 
assessing the merits of the theory. 

The present work is a continuation of work described in a previous paper, Part I, 
on bubble formation in a viscous liquid. 

Numerous investigators, [1,4-7,9,12,13,15,16,23,28,30] have measured the volume 
of bubbles formed when air or other gas is blown into an inviscid liquid. These 
experiments fall into three classes: 


(i) With a very low gas rate, [1,5,6, 12,30], the mechanism of bubble formation 
is similar to the mechanism of drop formation in the ”drop weight” method 
of determining surface tension. These gas flow rates are much less than those 
considered in the present paper. 


(ii) When the bubbles were formed at the end of a tube a few millimetres in 


6. Introduction 


diameter, [4,9,28,30], or at an orifice connected to a small gas buffer vessel, 
the conditions approached those of constant gas flow during bubble formation. 


(iii) Where the gas buffer vessel was larger, [7, 13, 15, 16,23], the conditions were 
more nearly those of constant pressure, though with an intermediate sized 
vessel, the volume of the vessel is an important factor determining bubble 
size. 


No previous attempts appear to have been made to predict bubble volumes and 
frequencies from theory, though HAYES, HARDY, and HOLLAND, [13],have written 
down equations describing the vertical motion of the bubble during formation. 
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7. Experimental Methods 


7.1. Apparatus for bubble formation 


The apparatus has been described in detail in chapter 2. The liquid in which 
bubbles were formed was contained in a vertical tube of 14.7 cm internal diameter 
with liquid depths up to 15cm. This tube was mounted with its axis vertical, on 
top of a 45L drum which acted as a buffer, gas being fed steadily into the drum 
and passing out as bubbles through an orifice at the base of the cylinder containing 
the liquid. Orifices of three kinds were used, and are illustrated in Figure 2.2: 


1. For the constant pressure experiments, the orifice was essentially a plain hole 
in a horizontal plate with a diameter between 0.29 and 0.46 cm. 


2. For the constant flow experiments with large gas flow rates, the orifice 
consisted of two plates. The top plate was in contact with the liquid and was 
drilled with a plain hole with a diameter between 0.30 cm and 0.5cm. The 
bubbles formed immediately above this hole, and to ensure constant flow a 
sintered plate was brazed on below the hole. 


3. For the constant flow experiments with smaller flow rates, the gas was passed 
into the liquid through a long capillary. The end of the capillary where the 
bubbles formed was let into a horizontal plate so that the bubbles were formed 
under conditions similar to those with orifices (a) and (b) in Figure 2.2. 


The arrangements for measuring gas flow and drum pressure were described in 
chapter 2. 

In most of the experiments the bubble frequency was measured by a stroboscope. 
In some of the constant pressure experiments, the bubble formation was not regular 
enough for the use of a stroboscope, and it was necessary to take ciné photographs 
to measure the volumes of the individual bubbles. For this purpose a PATHE WEBO 
ciné camera was available, with up to 80 frames, and a minimum exposure time 
of about 1/600 s. 


8. Constant Flow Regime 


8.1. Theoretical Considerations 


The assumptions used in deriving the theory are similar to those described in the 
previous paper: the bubble is assumed to be spherical at all times during formation, 
and the upward motion is determined by a balance between the buoyancy force 
and the upward mass acceleration of the fluid surrounding the bubble. Thus the 
equation of upward motion is: 


d /11 ds 


V being the volume of the bubble at time t after it started growing, and s is the 
vertical distance of the centre of the bubble above the point where the gas enters 
the liquid. 

The assumptions and terms neglected in deriving Equation 8.1 are: 


(i) The flow around the bubble is assumed to be irrotational and unseparated, 
and the effective inertia of the surrounding fluid has been taken as HoV. 
Because of this assumption, the drag coefficient is zero. HAYES, HARDY, 
and HOLLAND, [13], assigned a drag coefficient to a forming bubble, thereby 
assuming that there is a fully established wake behind each bubble as it 
forms. However, when a solid sphere or cylinder is accelerated from rest 
in a fluid, the initial motion is practically irrotational, and the wake is not 
fully established until the body has moved an appreciable distance, [11,22]. 
It seems likely that the flow around the accelerating bubble will be similar, 
and consequently the drag coefficient is likely to be negligible until after the 
bubble has detached from the orifice. 


(ii) Neglection of the upward momentum of the gas leaving the orifice is justified 
by calculations described in chapter 3 


With a constant gas flow G, and an initial bubble volume of zero, V = Gt; by 
substituting this relation in Equation 8.1, the latter can be integrated. When 


t=0, ds = 0, and assuming that the bubble growth is terminated when its radius 


8. Constant Flow Regime 


r equals s the distance travelled, gives the final bubble volume ! 


2\ 5 

Vi = 1.378 - (E) l (8.2) 
g 

As each bubble breaks away from the orifice it leaves behind a volume of gas Vo 

which forms a nucleus for the next bubble. In the derivation of Equation 8.2 it was 

assumed that Vo was much less than V;, but with larger orifices Vo is not negligible 

and it is necessary to put 


4 
V = zur“ =Gt+V, (8.3) 
when integrating Equation 8.1. This gives, with the initial conditions s = ds = 0, 
169g (È Vot v Gt + Vo 
= : = =] . A 
a g E A m er 


The bubble is assumed to detach when s = r + ro, and the lifetime of the bubble is 
obtained by plotting Equation 8.4, and plotting r + ro as a function of time from 
Equation 8.3 on the same axes. Here ro is the radius of the orifice supplying the 


gas and it is assumed that Vo = ¿Tri 


8.2. Experimental Results 


The results are classified according to the flow rate; when this is small, Equation 8.2 
is appropriate, and Equation 8.3 and Equation 8.4 are used for larger flow rates. 
Figure 8.1 shows the results for small flow rates compared with Equation 8.2. It will 
be seen that above a flow rate of 1.5mLs”! the results are in excellent agreement 
with theory. At the lower flows the discrepancy between theory and experiment is 
greatest for the larger orifice with air bubbling into water. Approximate calculations 
suggested that the discrepancy is due to the influence of surface tension acting 
round the rim of the orifice, as it does in ” drop weight” experiments, and of course 
as the flow rate tends to zero the bubble volume is determined by the balance 
between surface tension and buoyancy forces, as in the "drop weight” method. 
Figure 8.2 shows the results for higher flow rates, compared with theoretical curves 
calculated from Equation 8.3 and Equation 8.4 using various values of rọ. These 
theoretical curves give bubble volumes that are somewhat higher than those given 
by Equation 8.2. With these higher flow rates, the double and quadruple bubble 
formation observed by HELSBy and Tuson, [14], was seen to occur. This made 
the stroboscopic method impracticable, and to measure the bubble volume it 
was necessary to take ciné pictures of the forming bubbles. The volume of each 


11.378 = 1. (22)* -112 
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Figure 8.1.: Experimental and theoretical bubble volumes for small constant flow 
rates 
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individual bubble could be calculated from the pictures, and for comparison with 
theory Figure 8.2 shows the volume of the leading bubble in the series of two or 
four, before coalescence. 

Figure 9.1 shows bubble formation in groups of two. It will be seen that the 

forming bubble is roughly spherical, and that each bubble does not assume the 
spherical cap shape, characteristic of fully separated flow, [8], until it has moved 
some distance from the orifice. These photographs therefore justify, to some extent, 
the assumption of unseparated flow during formation. 
Figure 8.2 shows good agreement between theory and experiment for flows between 
15 and 20mLs”!, as in the range 1.5mLs”!* to 3mLs”* shown in Figure 8.1. Data 
were not taken for the intermediate range 3.0mLs”* to 15mLs”! but there is 
good reason to suppose that theory and experiment would agree in this range. 
At flow rates above 20mLs”! the divergence between theory and experiment is 
believed to be partly due to the upward current induced by the bubbles in the 
liquid surrounding the orifice and partly due to the deformation of the forming 
bubbles from the spherical shape. Because of the upward current, each forming 
bubble is dragged upwards and therefore has a smaller volume than the theoretical 
bubble which is imagined to form in a stagnant liquid. Also, the flattening of the 
base of the forming bubbles may cause them to detach earlier than the idealized 
spherical bubbles. 
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Figure 8.2.: Experimental and theoretical bubble volumes for large constant flow 
rates 
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9. Constant Pressure Regime 


9.1. Theoretical Considerations 


The upward motion of the forming bubble is determined by Equation 8.1 in the 
same way as for the theory of the constant flow experiment. As described in 
subsection 3.2.1, it is supposed that gas flows into the bubble, from a vessel at 
constant pressure, through an orifice whose characteristic constant k is determined 
by an experiment in which the gas under consideration flows steadily through the 
orifice in the absence of liquid. The orifice is such that the volume flow rate of gas 
is k-times the square root of the pressure difference across the orifice, and when 
the bubble is forming, the flow rate is given by 


dV 20 
AP Ba ; 
J + pgs a (9.1) 


Here P = P, — pgh. Equation 8.1 and Equation 9.1 were solved numerically by 
the computer Edsac 2, with the initial conditions at t = 0: 


S == 
ds = 
dt 
4 
Vo 370 


The bubble was assumed to detach when s = r + rọ, and the results are plotted in 

the form of dimensionless groups in Figure 9.2, Figure 9.3, and Figure 9.4. The 

dimensionless groups were formulated from Equation 8.1 and Equation 9.1 with 

o and ro = 0, when there are unique relations between the dimensionless bubble 

volume V’ = -“,, the dimensionless pressure P’ = —>—, and the dimensionless 
pik2 gptk? 


flow G’ = el Here V = V, — Vo, and G is the mean flow rate during bubble 
92 p 


formation, and is V/t,, where t; is the time at detachment. For the case when 
a =0 and rọ = 0, these variables are related by the curves shown in Figure 9.2, 
Figure 9.3, and Figure 9.4, and the diagrams also show results for various values of 
o and ro. It is apparent that the effect of o and ro upon bubble volume and mean 
gas flow rate is small. Surface tension does however affect the minimum value of P 


for bubbling to occur, since bubbling will cease if P is less than me It should also 
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Figure 9.1.: Bubble formation in groups of two with a constant flow rate of 
13.7mLs!. ro =0.2cm 
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be noted that of Figure 9.2, Figure 9.3, and Figure 9.4, only two are independent, 
and the third can be plotted from the other two. In most practical problems the 
gas flow rate is fixed beforehand, and the bubble volume V and the pressure P 
have to be determined. 

The significance of the diagrams can be better understood by considering two 
limiting cases as follows: 


9.1.1. k SMALL; APPROACH TO CONSTANT FLOW 


When the orifice constant k is very small, P will become large for any flow to pass, 
and the flow will not vary during bubble formation. Substituting G = k - VP in 
Equation 8.2 it is easily shown that V’ = 1.378P"5, G’ = VP’, and these relations 
are shown plotted on Figure 9.2, Figure 9.3, and Figure 9.4, and, as would be 
expected, they approach the ” constant pressure” curves asymptotically. 

In Figure 9.3, the horizontal distance between the ” constant pressure” curve and 
the "constant flow” line is a measure of the difference between the actual pressure 
drop through the system and the pressure drop which would be obtained if the 
flow through the orifice were uniform and equal to G. This means that the normal 
method of calculating pressure drops through bubbling devices by adding together 
a static liquid head, a surface tension term, and a hole pressure-drop term gives 
a total pressure drop which is too high, since all these pressure drops do not act 
simultaneously. 
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Figure 9.2.: The theoretical relations between dimensionless bubble volume and 
pressure. When P’ = 0, V’ = 0.500 for constant pressure. The points as well as 
the curves were derived from theory 
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Figure 9.3.: The theoretical relations between dimensionless flow rate and pressure. 
When P’ = 0, G' = 0.351 for constant pressure. The points as well as the curves 
were derived from theory 
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Figure 9.4.: The theoretical relations between dimensionless bubble volume and 
flow rate. At the point X, V’ = 0.500, G’ = 0.351 
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9. Constant Pressure Regime 


9.1.2. P=0; THE CRITICAL FLOW 


The ” constant pressure” curves in Figure 9.2 and Figure 9.3 show that when P = 0, 

G and V have finite values. This discontinuity is due to the presence in Equation 9.1 

of the term pgs, which gives a finite pressure difference across the orifice as soon as 

the bubble starts to move. With a finite value of the surface tension, the minimum 
20 


value of P is a) and if V is plotted as a function of P, the curve will show a step 
up, at P = os from zero to the value of V given by Figure 9.2. 


With o = 0 and P = 0, Equation 9.1 becomes 


dV 
Voss, 


and this equation together with Equation 8.1 has the following analytical solution 
for rg = 0: 


2p 2 
Zul, 9.2 
V 33 kgt (9.2a) 
8 2 
= — gt”. 9.2b 
¿230 (9.2b) 


The bubble will detach when s = (295, giving the volume V of the detached 
bubble and the mean flow rate G as follows: 


y 3 1 
V 33\4 (3 \2 
Ve Zee Se) =0.500 9.3 
k3 pi (=) (=) i 0) 
= 1 1 
2 8 (3\4 
N -L (3) (2) = 0.351. (9.4) 
g2k1p8 4 2 At 


These numbers agree well with the values calculated from the computer results. 


9.2. Experimental Results 


Experimental results with the "constant pressure” arrangement are shown in Ta- 
ble 9.1, and Figure 9.5, Figure 9.6, and Figure 9.7 show the results graphically for 
To = 0.149cm. These results show clearly the phenomenon of the critical pressure, 
with a sudden rise in bubble volume and mean flow rate when P is nearly m 


Comparing Figure 9.5 and Figure 9.6, it can be seen that when G is reduced below 
its critical value (about 32mLs~‘), V changes very little, and this is in accordance 
with the theory which predicts that V should stay constant when G is less than 
its critical value. For this reason, the theoretical line on Figure 9.5 is shown chain 
dotted, since when P = 968, V can be either zero or 3.5. But in Figure 9.6, the 
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Figure 9.5.: Bubble volume as a function of pressure with air and water. k = 1.90, 
ro = 0.149 cm 
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Figure 9.6.: Mean flow rate as a function of pressure, with air and water. k = 1.90, 
ro = 0.149 cm 


theoretical line is shown continuous, because at P = 968, all values of G between 
0 and 67 are possible. Intermediate values of G between 0 and 67 give a bubble 
volume of 3.5mL. Figure 9.5 and Figure 9.6 show that when P is greater than 
968 the theoretical flow rate and bubble volume are considerably greater than the 
experimental values, and the differences are thought to be due partly to the upward 
current of water caused by the stream of rising bubbles, and partly due to bubble 
deformation; these two factors were discussed in connection with constant flow 
results shown in Figure 8.2. 

Figure 9.7 shows bubble volume V plotted as a function of mean flow G. When 
G is less than 67, the theoretical value of V is constant at 3.5, and the dotted 
curves are merely continuations of the theoretical curves for ” constant pressure” 
and ”constant volume”. At a given value of P, the experimental values of V and 
G are both less than the theoretical, and consequently the points on Figure 9.7 are 


59 


9. Constant Pressure Regime 


Theory 


Critical bubble volume 


SS 
S 0:0222G* 


V - BUBBLE VOLUME (m1) 


0 20 40 60 80 
G- MEAN FLOW (mi/s) 


Figure 9.7.: Bubble volume as a function of mean flow for constant pressure with 
air and water. k = 1.90, ro = 0.149cm. X is the experimental point of critical 
pressure 


56 


9. Constant Pressure Regime 


0.187cm 


-l 
2g 2, Tp = 


1 


26 mLs™t, k = 


G= 


9.8.: Bubble formation at constant pressure with a mean flow rate below 
3.06 mL cm 
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9. Constant Pressure Regime 


Figure 9.9.: Bubble formation at constant pressure with a mean flow rate above 
the critical value. G = 60mLs”!, k = 3.06 mLcm2 g”2, ro = 0.187 cm 


shifted down towards the bottom left hand corner of the diagram, and this shift 
includes the critical flow point X. In design work, G is the independent variable, 
and Figure 9.7 shows that for low values of G, the actual value of V is less than 
the theoretical, but for high values of G the simple formula from the ”constant 
volume” theory gives an approximate representation of the actual bubble volumes. 
Table 9.1 shows that the characteristics of Figure 9.5, Figure 9.6, and Figure 9.7 
are accentuated for the larger orifices. 


98 


Mean flow G 


Bubble volume V 


a HI, 

k ro 20 P Expt. Theory Expt. Theory = 
(mL em? ga) (cm) (g Pe s2) (gcm (mL s7!) (mL s”?) (mL) (mL) IDEE et 
1.90 0.149 968 951 17 67 2.0 3.5 0.67 
1.90 0.149 968 951 23 67 2.1 3.5 0.96 

1 0.90 0.149 968 951 32 67 2.3 3.5 1.42 
1.90 0.149 968 1118 45 70 2.9 3.8 2.14 
1.90 0.149 968 1323 61 76 3.4 4.2 3.06 
3.06 0.187 771 779 26 102 3.2 6.1 1.11 
3.06 0.187 771 779 33 102 3.4 6.1 1.48 
3.06 0.187 771 877 47 105 4.1 6.4 2.26 
3.06 0.187 771 1024 60 112 4.5 6.9 3.04 
3.82 0.206 698 734 17 124 4.1 7.8 0.67 
3.82 0.206 698 734 30 124 4.3 7.8 1.31 
3.82 0.206 698 832 57 129 4.9 8.3 2.84 
3.82 0.206 698 1006 68 141 5.7 9.1 3.53 
4.90 0.230 625 632 25 156 5.6 10.7 1.06 
4.90 0.230 625 739 60 163 6.9 11.4 3.04 
4.90 0.230 625 790 68 169 7.1 11.7 3.53 
4.90 0.230 625 800 70 169 7.5 11.8 3.64 


Table 9.1.: Formation of Air Bubbles in Water with Constant Gas Pressure o = 72 dyn cmt 


9. Constant Pressure Regime 


In comparing ” constant flow” and ” constant pressure” experiments (Figure 8.2 
and Figure 9.7), it should be noted that in the latter experiments the effect of the 
”neck” formed during detachment will be more important, because in the one case 
V is proportional to t and in the other case to t?. This means that the bubble 
formed in the "constant pressure” experiment will receive a greater proportion of 
its gas during the detachment period. Figure 9.8 and Figure 9.9, which are ciné 
pictures of bubbles forming in a "constant pressure” experiment, do show the very 
great rate of bubble expansion during detachment, compared with the ” constant 
flow” bubbles shown in Figure 9.1. Figure 9.8 and Figure 9.9 also demonstrate 
bubble formation with pressures above and below the critical value. Figure 9.8 
shows that, below the critical pressure, there is a discrete time interval between 
one bubble and the next — the magnitude of the interval depending upon the size 
of the gas reservoir below the orifice — whereas Figure 9.9 shows continuous bubble 
formation. 
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10. Conclusions 


Experiments with a constant gas flow rate during bubble formation have shown 
that Equation 8.2 can be used as a rough guide to the relation between the bubble 
volume V and the gas flow rate G. This equation can also be used when the bubbles 
form above an orifice with a constant supply pressure below, provided the gas flow 
rate is somewhat above the critical flow rate. 

If the gas flow rate is less than or equal to this critical flow rate, the bubble size 
will be roughly constant, and the time interval between the bubbles will depend 
upon the volume of the vessel below the orifice. If this vessel is infinitely large, 
bubble formation will cease altogether when the pressure P is less than — 

It seems likely that the critical flow rate described in the present paper is closely 
related to the flow rate at which downward leakage or ”dumping” will occur on a 
sieve plate, because if the total gas flow divided by the number of holes is less than 
the critical value of G, some of the holes will not be bubbling continuously, and 
the opportunity for leakage will occur. 
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A. Effect of the Perspex ring on 
bubble frequency 


At the outset of the investigation it was decided that the bubbles should be formed 
in a as stagnant as possible region of the liquid as possible. It is clear that if 
bubbles are periodically formed at a rapid rate, the liquid will circulate and this 
will influence the bubble volume. This effect is difficult to account for theoretically, 
and a ring to reduce circulation, in the region where the bubble is actually formed, 
was therefore fitted to the orifice as shown in Figure 2.2. The effect of such a ring 
is best illustrated by means of an example. 

Air at the rate of 12.5 mLs”* was passed into a liquid of viscosity 480 cP (480 mPas), 
and the frequency of bubble formation determined for various diameters of the 
Perspex ring. Figure A.1 shows a graph of bubble frequency versus ring diameter 
from which it can be seen that the frequency approaches an almost constant value 
when the ring diameter is greater than about 1.5 times the bubble diameter. 

All the results quoted in the present paper were obtained using rings with diameters 
at least one and a half times as large as the diameters of the bubbles. In each case 
a preliminary experiment was carried out to ensure that the results conformed to 
the flat section of the curve. 


A. Effect of the Perspex ring on bubble frequency 
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Figure A.1.: Bubble frequency as a function of orifice ring diameter 


Gas flow rate = 12.5mLs”! 

Liquid viscosity= 480 cP (480 mPas) 
Orifice radius = 0.0975 cm 

Liquid seal = 7.9 cm 


65 


B. The effect of the gas momentum 
on bubble size 


In the experiments described above, gas was blown from a small orifice and it thus 
had a certain amount of upward momentum which was neglected in the theoretical 
treatment. The force due to the rate of change of momentum is en. where Ap 
is the effective area of the orifice and pa the density of the gas. In the present 
experiments this force acts in the same direction as the buoyancy force, and, if it 
is to be taken into account, the equation of motion of a bubble forming at small 


constant flow rates can be rewritten in the form: 


(B.1) 


i 
Rearranging and substituting r = ($2) 3 gives 


TT 


ds 29 (3Gt\? Ga (4rY? 
dt 9v \ 4r = 6rAopv \3Gt) ” 


which upon integration and substitution of V/G for t gives 


¿2 3GY5 (V y G?pa An\3 (UN? (B.2) 
iby \ 4a G)  AnAopv \3G GJ ` l 
If the bubble is assumed to detach when s = r then from Equation B.2 the condition 
prevailing at detachment will be 


1 2 
29 3 \3 4 Gpe An \ 3 1 
lis -(—] -V34 . V3. B. 
15vG (2) us Ar Agpv ( 3 ) (Big) 


The first of the two terms on the right hand side of Equation B.3 is due to viscous 
forces and the second to momentum forces. If G = 1.0mLs”?, v = 640 cS (u = 
800cP), p = 1.25gmL~", pg = 0.0012g mL! and Ay = 0.0035cm?, then the 
momentum accounts for only 0.5% of the right hand side and its effect can thus 
be neglected. 

For high gas flow rates Equation B.1 has to be modified to allow for the effect of 
liquid inertia. If G = 17mLs”*, v = 586cS (u = 735cP), p = 1.2548 mL ?, and 


win 


B. The effect of the gas momentum on bubble size 
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Figure B.1.: Bubble sizes compared for carbon dioxide and air 


Flow rate = 17 mL s7! 
Liquid viscosity = 735 cP (735 mPas) 
Orifice radius = 0.0745 cm 


Ay = 0.0175 cm?, then the theory predicts that the volume of an air bubble is 1.0% 
greater than a carbon dioxide bubble. The theory has been tested experimentally for 
this viscosity and Figure B.1 shows that for a flow rate of 17mLs”! the air bubbles 
are 1.8% larger than the carbon dioxide bubbles. In the theoretical calculations, 
the actual and effective orifice areas were assumed equal. For constant flow systems, 
the difference in bubble size can thus be accounted for in terms of the momentum 
of the gas. For constant pressure systems the effect of the density of the gas on the 
orifice discharge coefficient must also be taken into account. 
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C. Pressure necessary to generate a 
bubble in a viscous liquid 


Suppose that a gas is being blown from a source into an infinite liquid in the 
absence of gravity. The resulting bubble will be spherical, and its centre will be 
stationary. Under these conditions PORITSKY, [21], has shown that the viscous 
terms vanish from the NAVIER-STOKES equation, and therefore the mean pressure 
within the liquid at any point is the same as it would be in an inviscid liquid. 
However, the pressure within a viscous liquid in motion is not the same in all 
directions, and at the surface of the bubble the radial pressure differs from the 
mean pressure by ur : z [21]. This term gives the amount of the increase of 
pressure within the bubble due to the viscosity of the liquid, and in the constant 
pressure experiments the flow into the bubble would be reduced by including the 
term within the square root in Equation 3.12. 

The last column in Table 3.4 gives the term: 


4pv dV 20 
=1-,/1- . dl ae el 
(Pre) 


which is the fractional variation in the flow into the bubble caused by liquid viscosity. 
These figures show that although the liquid viscosity has a very large effect at the 
start of the motion, the overall effect on the flow is likely to be about 6%. From 
Table 3.3 it follows that the effect on the final bubble volume is likely to be about 
3%. 


D. Pressure due to kinetic energy 
imparted to the liquid 


If a gas is blown into a liquid the latter will acquire kinetic energy due to the 
expansion of the resulting bubble. The total kinetic energy generated during the 
formation of a bubble is given by the equation: 


Kinetic energy = f 2pru*RidRs . (D.1) 


Substituting E - R? for u and integrating we have 


2 
Kinetic energy = Ce ; (D.2) 
nr 
P, the mean pressure due to the kinetic energy of the radial motion, is obtained by 
equating the work done during bubble formation PGt to the kinetic energy given 
by Equation D.2. Calculations showed that, for the cases considered in the present 
paper, P is negligible in comparison with the other pressures involved. 


E. Note on the formation of bubbles 
in an inviscid liquid 


In the discussion on the formation of bubbles at small constant flow rates in liquids 
of high viscosity it was assumed that inertia forces were negligible in comparison 
with viscous forces. For the case of an inviscid liquid, the viscous forces are 
negligible in comparison with the inertia forces. 

The equation of motion is then 


d (11 ds 
Mo= EM E.1 
I~ at E =) a 
which on integrating twice, using the condition that at t = 0, ds = s = Q, gives 
4g 15 
=—.-f. E.2 
a (E.2) 
1 
At detachment r = s = (2+) 3 and hence the final bubble volume V is given by 
zye /11\8 GÊ 
V=Gt=(—) -(—) -G3 = 1.378. 2 (E.3) 
Ar 4g g5 


It is interesting to note that, except for the numerical constant, Equation E.3 is 
identical with the empirical expression given by VAN KREVELEN and HOFTIJZER, 
[30], to predict the volume of a bubble in terms of gas flow rate. VAN KREVELEN 
and HOFTIJZER found that, provided the gas density was negligible in comparison 
with the liquid density, the volume of a bubble was given by the expression: 


725 Gs Gs 
=. =172%. 


V = 
675 g5 g5 


(E.4) 


F. John Frank Davidson. Obituary 
Notes 


JOHN FRANK DAVIDSON, FRS, FREng, FIChemE, (7 February 1926 — 25 December 
2019) was a British chemical engineer and former Shell Professor of Chemical 
Engineering at the University of Cambridge. He is regarded as the founding father 
of the subject of Fluidization in Chemical Engineering. 

DAVIDSON’s works made an important contribution to the development of the 
hydrodynamics of a circulating fluidised bed and the heat transfer in a fluidised 
bed and also to the creation and implementation of methods for lignite combustion 
in a fluidised bed. DAVIDSON received the degree of PhD in 1953 and the degree 
of Doctor of Science in 1968 at the University of Cambridge. In 1974, for works on 
two-phase flows and, first of all, for achievements in fluidisation, he was elected 
a Fellow of the Royal Society in chemical engineering (He was elected as a Vice- 
President, Royal Society of London, 1989). In 1974-1975, he was a member of the 
Court of Enquiry for the Flixborough disaster. 

DAVIDSON’s research work was always closely related to his pedagogical activities. 
He started as a university demonstrator in chemical engineering (1950) and then 
was a university lecturer (1954), a reader in chemical engineering (1964), and a 
professor (1975). In 1975--1993, DAVIDSON headed the Department of Chemical 
Engineering of the University of Cambridge and did much for the development of 
the department, updating of the curriculum, and strengthening of relations with 
industry. In 1978-1993, he bore the title of Shell Professor of Chemical Engineering 
since Shell was an important sponsor of the foundation of the department. 
Indissoluble are the ties of DAVIDSON’s life and activities with Trinity College, one 
of the oldest Cambridge colleges, ties lasting from 1949 to the present time. In 1957, 
he became the College’s steward responsible for the entire household, including 
organisation of receptions of the royal family. The most difficult challenge in this 
period was the reconstruction of the Old Kitchen, built in 1605. This reconstruction 
lasted several years in the early 1960s. In 1992-1996, DAVIDSON was vice master 
of Trinity College. 

In 1970-1971, DAVIDSON was president of the Institution of Chemical Engineers; 
he is a Foreign Associate of the US National Academy of Engineering (1976), 
Docteur Honoris Causa de l'Institut National Polytechnique de Toulouse (1979), an 
Honorary Doctor of Science of the University of Aston (1989), and a Foreign Fellow 
of the Indian National Science Academy (1991). He has the Kurnakov Memorial 
Medal given by the Kurnakov Institute of General and Inorganic Chemistry of 


F. John Frank Davidson. Obituary Notes 


Figure F.1.: JOHN FRANK DAVIDSON, bromide print, 1993 


the USSR Academy of Sciences (1991) and was elected a foreign member of the 
Russian Engineering Academy (1998). In 1999, DAVIDSON was awarded the Royal 
Medal of the Royal Society of London. 

DAVIDSON was a frequent participant at and guest of MENDELEEV congresses, at 
which he has delivered plenary lectures. He continued active research work at the 
Department of Chemical Engineering of Cambridge University after his retirement 
in 1992. 

He was awarded a Royal Medal in 1999 for his chemical engineering work, and has 
been a Fellow of the Royal Society since 1974. He was president of the Institution 
of Chemical Engineers (IChemE) from 1970 to 1971 and vice-president of the Royal 
Society in 1989. He was appointed a Fellow of the Royal Academy of Engineering. 
In 1994, the National Portrait Gallery, London, along with the Royal Society, 
commissioned a portrait of DAVIDSON by photographer NICK SINCLAIR, which 
was acquired for its permanent collection. 
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